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ABSTRACT 

The in-line and transverse forces acting on cylinders of 

1.0 inch to 2.5 inches diameter placed near a plane wall in a 

harmonically oscillating flow have been measured. The drag 

and inertia coefficients C, and C for the in-line force and 

d m 

the lift coefficients and C^ T for the transverse forces 

away and toward the wall, respectively, have been determined 
as a function of the relative gap e/D between the wall and the 
cylinder and the period parameter V^T/D. The relative gap 
ranged from 0.0143 to unity and the period parameter from zero 
to 35. In the subcritical range where these experiments have 
been performed, the effect of the Reynolds number, which ranged 
from 4,000 to 30,000, was found to be secondary and certainly 
obscured by the excellent correlation provided by the period 
parameter . 

The frequency of the transverse force was found to range 
from one to five times the frequency of oscillation of the har- 
monic motion, but a dominant frequency of twice the oscilla- 
tion frequency was seen over most of the values of V^T/D. 
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NOMENCLATURE 



amplitude of the motion 
average drag coefficient 

maximum lift coefficient away from the boundary 
maximum lift coefficient toward the boundary 
average inertia coefficient 
diameter of test cylinder 

normal distance between the boundary and the cylinder 
surface 

instantaneous total force acting on the test cylinder 

drag force acting on the test cylinder 

lift force acting on the test cylinder 

water depth 

wave number (2 tt/1) 

length of test cylinder 

wave length 

period of oscillation 

time 

horizontal fluid particle velocity 
instantaneous velocity 
maximum velocity 

vertical fluid particle velocity 

vertical position of the test cylinder 

percent error 

fluid kinematic viscosity 

fluid density 

wave frequency (2 tt/T) 
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I. 



INTRODUCTION 



A. EXPERIMENTAL JUSTIFICATION 

With the current energy deficit becoming of vital impor- 
tance, offshore oil production assumed great priority, and 
pipe lines became the most feasible mode of liquid transporta- 
tion to the shore. This, in turn, necessitated the deter- 
mination of the wave and/or current induced forces on submerged 
pipes placed at or near the ocean bottom. The problem, how- 
ever, is not unique to the pipe lines. 

Underwater cables for communication and sensing purposes 
and all submerged structures, be they offshore oil platforms 
or oceanographic research instruments, are subjected to wave 
induced forces . 

The proximity of pipe to a solid boundary, such as the 
ocean bottom, further complicates the problem, and, unless 
the submerged pipe lines are anchored to a buoy system or 
buried, the bottom effects must also be considered. Even 
buried pipe lines, subject to scouring can become uncovered 
and have their support dug out from under, leaving them sen- 
sitive to wave induced forces. 

B. BACKGROUND THEORY 

The in-line force acting on a cylinder immersed in a time- 
invariant flow may be expressed, according to Morison and his 
co-workers [6] as: 
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( 1 ) 



F=F ,+F = 1/2LC ,DpV 2 +C ttD 2 /4 LdV/dT 

dm a m 

in which and represent time-invariant drag and inertia 
coefficients, respectively. The first term on the right-hand 
side of this equation represents the velocity-squared depend- 
ent drag force and the second term the acceleration-dependent 
inertia force. 

A simple dimensional analysis shows that the time- 
averaged drag and inertia coefficients for a cylinder immersed 
in a harmonically oscillating uniform flow, represented by 

V=-V cos 2irt/T 
m 

depend on the period parameter V T/D and the Reynolds number 
V m D/v. In addition to the in-line force, the cylinder is 
subjected to an alternating transverse force due to separa- 
tion and vortex shedding. The proximity of a solid boundary 
can alter significantly the magnitudes of both the in-line and 
transverse forces and requires the determination of the wall 
proximity effect on the drag, inertia, and lift coefficients. 
In fact, it is the determination of this effect that prompted 
the present study. 

A harmonically oscillating uniform flow represents only 
approximately the oscillatory motion induced by waves about 
a submerged cylinder. Nevertheless, the wave-induced forces 
may be predicted with sufficient accuracy through the use of 
drag, inertia, and the lift coefficients obtained with har- 
monic flow, provided that the velocities and accelerations 
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about the cylinder are either measured or calculated through 
the use of an appropriate wave theory. 

For small amplitude waves, Airy's theory may be used to 
evaluate the velocities and accelerations of a given time and 
depth. The limitations of this theory are well-known and will 
not be represented here. 

The horizontal and vertical components of velocity may be 
written as [11] , 

u=Ao(cosh k(h+z)/sinh kh) cos 0 (2) 

w=Ao(sinh k(h+z)/sinh kh) sin 0 (3) 

in which A represents the wave amplitude, k the wave number 
defined by 2n/l, h the water depth, z the vertical position 

2 TT 

of the cylinder, a the wave frequency given by and 

0=kx -at . 

The components of acceleration are given by 

6u/<5t=-Aa^ (cosh k(h+z)/sinh kh) sin 0 (4) 

2 

6w/ St=Ao (sinh k(h+z)/sinh kh) cos 0 (5) 

The in-line force may be expressed through the combination of 
equations (1) , (2) , and (4) to yield, 

2 2 

F=Aa LC D cosh k(h+z) sin 0 

m 4 sinh kh 

2 2 

-A C^LpD (cosh k(h-t-z) | cosat | cosat (6) 

(sinh kh) ^ 

The transverse force cannot, however, be calculated in a 
similar manner through the use of equations (1), (3), and (5). 
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It is fundamentally the eddy behavior that determines the 
transverse force. Thus, the lift resulting from a nonlinear 
motion is not susceptible to estimation by superposition of 
individual responses to each harmonic acting in isolation. 

The lift coefficients in the present study are evaluated 
as follows: 

C^ A = ( maximum transverse force away from the wall in a cycle ) 

0.5pDLV m 2 (7) 

C^ T = ( maximum transverse force toward the wall in a cycle ) 

0.5pDLV m 2 (8) 

Force coefficients based on other time-invariant magnitudes 
such as the root-square values could have been easily obtained. 
It will suffice here to remark that a simple Fourier analysis 
for the in-line force and the peak value for the transverse 
force proved adequate for all purposes. 

The determination of the coefficients con- 
stitute the essence of the present study for a cylinder placed 

within the proximity of a plane wall in a harmonically oscil- 
lating flow. 

C. RELATED RESEARCH 

Yamamoto, et al . , in 1974, studied wave forces on cylinders 
near a plane boundary in the range of Reynolds numbers from 
2,000 to 30,000 and the surface period parameter from about 
0.3 to 3. The total force was considered as the sum of 
components due to water particle velocity squared (lift and 
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drag forces) , and due to water particle acceleration (inertia 
force) . They have correlated the force coefficients and 
with the relative distance of the cylinder from the bound- 
ary. 

The vertical and horizontal hydrodynamic forces were 
obtained for unseparated flow over the submerged cylinder by 
using the method of double images and applying Blasius' theorum. 

The results indicated that the proximity of a plane 
boundary modified both the lift and the inertia coef f icients . 
However, these results are applicable only to unseparated 
flows where drag is negligible compared to inertial forces. 
Within a narrow range and relatively low values of the period 
parameter, it was shown that the maximum horizontal force oc- 
curred at zero crossings of the surface wave, indicating 
that wake dependent drag force is negligible and the horizon- 
tal force is composed mainly of the inertia force. The verti- 
cal force, for larger values of the relative distance was 
unsinusoidal with large downward forces (toward boundary) at 
the surface wave crests and small upward forces at surface 
troughs . 

In 1967, a theory for simple shear flow across a circular 
cylinder in proximity to a plane wall, was developed by Arie 
and Kiya [2]. Theoretically, they showed that a force acts 
downward (negative uplift) for uniform velocity, but a force 
acting upwards (away from the plane) exists for flow having a 
velocity gradient. Their experimental data failed to enforce 
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this theory, but did show that lift on the cylinder increased 
as the clearance between the cylinder and wall decreased. A 
wind tunnel was used for their experimentation. 

Studies on the effects of ocean currents on pipes anchored 
just above the ocean floor were conducted by Wilson and Cald- 
well [12] in 1970, and the lift and drag coefficients were 
determined. The frequency of the forces on the pipe was 
found to be effected by the eddy-shedding frequencies. 

In 1971, a study by Grace [3] on the effect of the clear- 
ance of pipe above a flume and the effect of orientation of 
the pipe to the wave crest was made. With clearance ranges 
from 1/32 inch to 1-1/2 inches and wave periods of 2 to 6 
seconds in 3 foot deep water, he found that the horizontal 
force normal to the pipe was not effected by bottom clearance 
and that this force decreased as the pipe became perpendicular 
to the wave fronts. He also found that the "vertical force" 
decreased as the bottom clearance decreased. 

In much of the work conducted on wave forces, correlations 
between C m , C^, and the Reynolds number V m D/v, where v is 
the kinematic viscosity, have been sought with little success. 
Wiegel, Beebe, and Moon (1957) [11] found no relation between 

these parameters in studies of ocean wave forces on cylindri- 
cal piles. Jen [4], in 1968, also failed to correlate C or 
with Reynolds number in a laboratory study of a 6 inch 
diameter pile. 

In 1958, Keulegan and Carpenter [5], investigating the 
inertia and drag coefficients on circular cylinders, found 
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a correlation between these force coefficients and a dimen- 
sionless parameter, V^T/D. Working in a rectangular basin 
with standing waves, they covered V^T/D values between 2.7 and 
120 over a Reynolds number range from 4000 to 29,300. The 
lowest value of and the highest value of were found to 

occur at a V T/D of about 15. 
m 

Sarpkaya and Tuter [8], in 1974, with an oscillating U- 
shaped channel and various cylinders, reconfirmed Keulegan- 
Carpenter data except for a small spread in the upper regions 
of V^T/D, where the Sarpkaya-Tuter data may be more accurate 
and reliable. An important finding of the Sarpkaya-Tuter work 
was that the lift or transverse force was as larage as the in- 
line force. 

The close correlation between C and C, and the period 
parameter V^T/D is evidenced by the close fitting curves re- 
lating these parameters. The curve of Kuelegan-Carpenter 
forms a lower envelope for ocean data gathered by Wiegal , Beebe, 
and Moon, while the curve forms an upper limit for the same 
data [11]. 

The work on wave induced forces on cylinders has attracted 
much attention, but results seem quite scattered and some- 
times contradictory. Also, the effect of a plane boundary in 
proximity to the cylinder presents a new variable. It is the 
purpose of this thesis to analyze the forces on a cylinder in 
the proximity of a plane wall and to correlate the various 
coefficients with the period parameter V m T/D, and the relative 
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distance between the cylinder and the wall, through the 
of harmonic motion. 



use 
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II. METHOD OF ANALYSIS 



To analyze the time dependent force in unsteady periodic 
flow, a Fourier analysis will be made, beginning with the 
Morison's equation given by 

F=l/2 C d pAV|v|+ C m pV dV/dt (9) 

2 

The enclosure of one of the components of V term in "absolute 
value" brackets accounts for the change in sign of the resist- 
ing force as velocity changes direction in an oscillating 
flow. 

For a circular cylinder, the equation reduces to: 

F=C ir/4LD 2 pdV/dt + ^dLDpV|v| (10) 

m 2 

Representing the velocity of the harmonic motion by 

V=-V cosat (11) 

m 

and inserting into equation (10) , one has 



2 _ 

= ir D 



pLV D/2 V T 
m m 



C sinat - C, cosot cosat 
m d 1 1 



(12) 



Multiplying both sides by sin t and integrating, one has, 



j ' F singt d ( pt) 
pV m 2 D/2 . L 



2tt 2 

0 / ^ C_sinatsinat d(at) 



V T nf 
m 



2tt 



- 0 / C, | cosat | cosat | sinat d(at) (13) 
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which gives 



2n 2 

r F sinat d(ct) = tt D C m (TT) 

° pLV 2 D/2 V T 

K m in 



(14) 



and 




2 V t 
■ m — 
tt 3^ 



2ir 

of 



Fsinot d (ot) 
P LV m 2D 



(15) 



Multiplying both sides of the force equation by cosat, and 
integrating, results in 



2ii 

of — 

pLV 



F 

2 

m 



D/2 



cosat d(at) 



2tt 2 

= „/ it D C sinat cosat d(at) - 

V T m 
m 

(16) 
2 TT 

of | cosat | cosat cosat d(at) 



2tt 

of F 

P LV m 



2 



_cosat d ( at) 
D/2 



-8/3 C d 



(17) 



Finally, 

C, = -3/4 <3 (°t) (18) 

pLV D 
m 

For the analysis of these force coefficients a computer pro- 
gram was used (See Appendix A) . The determination of V m is 
made by considering the wave amplitude. A, and the period, T, 
resulting in 



V = 2ttA/T 
m 



( 19 ) 
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Thus the coefficients C and C, may be written as, 

m a 

C_ = 2T 2 E F sin (2l£) d(t/T) 

m n3 n 2 r . T 

D LAp 

C, = -3T 2 9 ZFcos d ( t/T) 

a 8pD7TLA 



( 20 ) 



( 21 ) 



As an alternate method of computing the force coefficients, 
a least-squares method was also employed. The method of least 
squares consists of the minimization of the error between the 
measured and calculated forces. Letting F represent the in- 
stantaneous measured force and F c the force calculated through 
the use of equation (10) , and writing 



2 2 
E = (F-F ) 

C 

2 2 

and dE /dC m =0 and dE /dC^=0, one has 

0 _ -8 2 r F | cosat | cosat d(at) 

°d ° J 2 

3u pDLV 

K m 



( 22 ) 



(23) 



and 



C = 
m 



2V T 




2 TT 

0 f F singt d(at) 

pV 2 LD 
m 



(24) 



It should be noted that the Fourier analysis and the method of 
least-squares yield identical values and that the values 
differ only slightly. 

The error between the measured and calculated forces, 
particularly in the neighborhood of the maximum forces, may 
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be further minimized by choosing the square of the measured 
force as the weighting factor in the least-square analysis. 
Thus writing 

E 2 =F 2 (F-F c ) 2 (25) 

2 2 

and dE /dC.=0 and dE /dC =0, one has 
' d ' m 



C „^5^3 ir j4^2— 

d pDLV 2 f 4 f l _f 3 f 3 
m 



(26) 



and 



T f f -f f 

C = — _5ii_3i2— 

in 3 TAn 2 f^ff-f^f, 
n LAD 41 33 



(27) 



The functions f^ are given by 

2n _ 4 2 tt ^ 

f.=/ F cos at d(at) , f = / F |cosat| cosat d(at) 

O ^ O 



2 

f^o-f F sinat cosat| cosat|d(at) 

2tt 2 2 2tt _ 

f^= 0 / F sin at d(at) , f^= 0 / F sinat d(ct) 



(28) 



Equations (27) and (28) may be shown to reduce to equations 
(23) and (24) by replacing F n in equations (28) by F n 2 and 
carrying out the necessary integrations in which F does not 
appear . 

Each wave cycle was divided into 36 time intervals of 
At=2. 86/36 seconds, (T=2.86). The force for each interval was 
found and incorporated into the computer program to calculate 
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the appropriate coefficients through the use of the three 
methods of analysis cited above. 
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III. EXPERIMENTAL EQUIPMENT 

The basic oscillating flow system consisted of a U-shaped 
vertical water channel with a channel cross section of 18 by 
20 inches shown in Figure 1. Oscillations were created by 
pneumatic pressure applied at a closed end of the channel. 

The closed end was opened by means of a slider crank mechanism 
uncovering a large exit and releasing the air pressure. The 
resulting oscillation was a near perfect harmonic as evidenced 
by the elevation and acceleration traces shown in Figure 2. 

The maximum amplitude of the oscillation was 11 inches and 
the natural damping of the oscillations was in the order of 
1/8 inch per cycle. The water level at .its minimum oscilla- 
tion height was twenty inches above the test cylinder. The 
cylinders were manufactured out of plexiglass tubes or rods at 
desired diameters and at lengths approximately 1/16 inch less 
than the channel width. Self-aligning bearings were imbedded 
at each end of the cylinders and caution was taken before each 
experiment to confirm the clearance of the cylinder end from 
the channel walls. 

The lateral and in-line (with respect to the fluid flow) 
force measuring devices consisted of various cantilever beams 
mounted to the ends of the test cylinder. Eight piezoresistive 
strain gauges were mounted on each cantilever beam and properly 
waterproofed. 

These force transducers were repeatedly calibrated by 
hanging known loads at the midsection of the cylinder in the 
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Fig. 2 Elevation 



and acceleration traces 
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horizontal and vertical directions. Some of the transducers 
were comprised of two pieces of thin cantilever beams cut 
orthogonal to each other and were able to measure both the 
in-line and lateral force simultaneously. Other transducers 
were capable of measuring only the in-line or lateral force. 
However, being rotatable -90°, this transducer could first 
measure the in-line force and with a 90° rotation, measure the 
lateral forces. In many cases the in-line and lateral forces 
were measured at each end of the cylinder and compared with 
each other. In no case did the force curves deviate from each 
other more than 5%, indicating a fairly uniform response along 
the cylinder with respect to the resultant forces. 

Throughout the investigation, the monitoring of the char- 
acteristics of the oscillations in the U-channel was of prime 
importance. Host of the difficulties in past determinations 
of and resulted from the difficulty of producing a 

purely harmonic motion or from having to determine oscilla- 
tion characteristics indirectly. The U-tube provides a per- 
fectly sinusoidal oscillation and the instantaneous displace- 
ment and acceleration are continuously recorded. The instan- 
taneous elevation in one leg of the channel was determined 
through the use of a capacitance wire connected to an amplifier- 
recorder system. Such wires have been used in the past to 
measure wave heights in open channels. The response of the 
wire was found through calibrations, to be perfectly linear 
within the range of oscillations encountered. 
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The instantaneous acceleration was measured by means of a 
differential-pressure transducer connected to two pressure 
taps at the mid-section of one of the channel walls. The 
acceleration was then calculated from 

Ap= psa (29) 

where A P is the differential pressure, p the fluid density, 
s the distance between pressure taps, and the instantane- 
ous acceleration. 

The effect of pressure drop due to the viscous forces was 
found to be negligible. 

The displacement and acceleration are, of course, in phase 
and may be used independently to calculate the velocity and 
displacement of the fluid, although displacement was actually 
used. The smoothness of the variation of acceleration and 
elevation traces shows the success in obtaining a purely 
harmonic motion as shown in Figure 2. 
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IV. DISCUSSION OF RESULTS 



The drag coefficients are shown in Figure 4 as a function 
of the period parameter V T/D, which also equals 2 ttA/D, for 
various values of e/D. The actual data points are shown only 
for one value of e/D, partly to simplify the presentation of 
the data and partly to give some idea of the variation of 
with respect to the three methods of analysis used in its 
evaluation. Firstly, it is apparent that all. three methods 
yield nearly identical results and that the data exhibit very 
little scatter. Secondly, the drag coefficient can reach 
values as high as 3.75. For e/D larger than unity, values 
approach those found for V m T/D=«. 

Evidently, it is not possible to explain the complex 
variations of with V T/D since it is largely determined by 
separation effects. 

The inertia coefficients are shown in Figure 5 in a manner 
similar to that for C^. The values of C rn corresponding to 
V m T/D=0 are obtained from the potential theory [10] . The ex- 
perimental values approach in all cases those predicted by 
the potential theory as V T/D->0. Once again, the inertia co- 
efficients approach those obtained for the limiting case of 
e/D=~ [ 8 ] for e/D values larger than unity. 

The lift coefficients (force away from the wall) are 

shown in Figure 7 and the lift coefficients (force towards 

the wall) are shown in Figure 6 . Also shown in Figure 6 are 
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Fig. 4 Drag coefficients versus period parameter 
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Fi g.7 Lift coefficients 



lift coefficients calculated through the use of the potential 
theory [10] for V m T/D=0. Aside from the fact that the lift 
coefficients can become very large for cylinders near a wall. 
Figs. 6 and 7 also show that the predictions of the potential 
theory for separated flow may be grossly in error. Figure 7 
shows that for a small gap (e/D=0 . 014 3) , large lift forces 
act on the cylinder in a direction away from the wall and C^ T 
is zero for all values of V^T/D larger than 15. According 
to the potential theory, a large net force exists toward the 
wall when a small gap exists between the cylinder and the 
wall [10]. Thus, the predictions of the potential theory can 
be used only for extremely small values of V m T/D for which 
separation does not occur. Evidently, the occurrence of separa- 
tion depends not only on e/D and V^T/D but also on the 
Reynolds number. No special effort has been made to determine 
the maximum values of V^T/D for a given set of e/D values 
and Reynolds number below which the predictions of the poten- 
tial theory may be applied. 

The practical significance of the occurrence of large 
lift forces both away from and toward the wall in a given 
cycle of oscillation is self evident. Such alternating forces 
could give rise to large vertical oscillations. These oscil- 
lations can cause scour and fatigue and lead to larger in- 
line forces than those calculated from equation (12) through 
the use of the coefficients from Figs. 4 and 5. In fact, it 
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is not difficult to imagine the occurrence of coupled in-line 
and transverse oscillations with disastrous consequences. 

The frequency f of the alternating transverse force was 
evaluated in terms of the frequency f of the fluid oscilla- 
tion and it was found, as done previously [8], that f is 
not a constant fraction of f and that several frequencies may 
occur during a given cycle. This seems to be rather reason- 
able because of the fact that in any cycle the magnitude of 
the velocity varies from zero to and that there can hardly 
be a single Strouhal frequency. In spite of the multiplicity 
of the frequencies of vortex shedding, the examination of the 
force records have shown that the vortices are shed with sig- 
nificant regularity over most of the V^T/D values at a domin- 
ant frequency of f v =2f. Other frequencies of the alternating 

lift force were f =f,3f,4f, and 5f. 

v 

Two additional matters need to be discussed, namely, the 
error between the measured and calculated in-line forces and 
the importance of the Reynolds number on the variation and 
correlation of the various force coefficients. 

The relative error in each cycle for each V T/D as well 
as the error corresponding only to the maximum in-line force 
was calculated for each method of data analysis. The detailed 
results which will not be reported here for the sake of brev- 
ity have shown that A defined by 
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X = (F (measured) -F (calculated) ) /F (calculated) 

max max ' max (30) 

is within -15% and reaches its maximum in the neighborhood 
of V m T/D=10 for the in-line force. This result may be taken 
as an indication of the fact that Morison's equation may be 
used to predict the in-line force acting on cylinders in the 
vicinity of a wall in a harmonically oscillating fluid through 
the use of the Fourier averaged coefficients. 

As to the Reynolds number defined by V^D/v, its role in 
wave force coefficients is not well understood. For Reynolds 
numbers below the critical value (a value which is expected to 
be considerably smaller than that corresponding to the steady 
flow about a cylinder, partly due to unsteady nature of the 
flow and partly due to the presence of the vortices both up- 
stream and downstream of the cylinder) , the Reynolds number 
appears to play some role in the variation of in the 
vicinity of the intermediate values of V T/D. For V T/D about 
10, only one or two vortices form on the downstream side of 
the cylinder during a given cycle. As previously observed by 
Sarpkaya [9], when the fluid accelerates rapidly at high 
Reynolds numbers, vorticity is slow to diffuse and therefore 
accumulated rapidly in the growing separation bubble behind 
the cylinder. Although this bubble of trapped vorticity 
soon reaches unstable proportions, the growth of the bubble 
and, hence, of the resulting vortices, are so rapid that the 



35 



vortices becoire much larger than their quasi-steady-state 
size. Furthermore, the relatively slow diffusion of vorticity 
(largely due to turbulence rather than molecular diffusion) 
enables the vortices to retain a larger fraction of the 
vorticity shed from the boundary layers as the vortices move 
a small distance downstream before the flow is reversed. Dur- 
ing this period, the drag on the cylinder rises above its 
steady-state value, both because the surface area covered by 
the separation bubble is temporarily greater than it would 
ultimately become in steady flow and because the enlarged 
wake distorts the outer flow, reducing the base pressure below 
its steady value. Thus, the effect of the Reynolds number, 
which is nearly absent in steady flow within the range of 
Reynolds numbers (4,000 to 30,000) discussed herein, is primari 
ly due to the differences in the rate of shedding and diffu- 
sion of vorticity between the steady and periodic flow. How- 
ever, the turbulent nature of diffusion and the excellent corre 

lation of C, with V T/D nearly obscure the effect of the 
dm 

Reynolds number, at least below the critical Reynolds number 
range. Near or above the critical Reynolds number range, how- 
ever, the effect of the Reynolds number may become more pro- 
nounced due to different reasons. It must be kept in mind 
that and are time averages and the error between the 
measured and calculated forces through the use of these time- 
invariant coefficients is not uniform and depends on V^T/D. 
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Similar arguments hold true for the transverse forces. The de- 
pendence of and C^ T on V^T/D is much stronger and tends 

to obscure the weaker dependence on the Reynolds number. 
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V. CONCLUSIONS 



Although the in-line and transverse forces acting on 
cylinders near a wall were measured for sinusoidally oscillat- 
ing flow rather than for wave-type fluid motion, the results 
presented herein enable a number of conclusions concerning 
the latter type of flow to be drawn. 

Morison's equation may be used to predict the in-line 
force acting on such cylinders or deeply submerged pipes (no 
free surface effects) subjected to local wave motion through 
the use of the Fourier averaged coefficients. The drag and 
inertia coefficients approach those predicted by the poten- 
tail theory as V^T/D approaches zero. 

The transverse forces, like the in-line forces, are 
strongly influenced by flow separation. Not only the magni- 
tude but also the direction of the forces differ significant- 
ly from those predicted by the potential theory. For intermed- 
iate values of the period parameter, potential theory under- 
estimates the transverse forces whereas for large values of 
V^T/D, it grossly overestimates them. 

Additional work with regard to the practical pipeline 
problems is needed to ascertain the effects of free surface, 
nonlinear waves, internal currents, roughness, and other 
environmental conditions . The data presented herein could 
form the basis for such extensions. 
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APPENDIX A 

(IN-LINE FORCE COEFFICIENT PROGRAM) 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

AVERAGE CD AMD CM CALCULATIONS FOR CYLINDERS 



TIME = DI MENSIONLESS TIME ( T I ME/ PE R 1 00 ) 

BETA = D! MENS I 3 ML: SS DISPLACEMENT ( UMAX * 3 E R/D I A ) 
AMP= A MP L I TUDE OF MOTION IN FEET 
Z I =F GRC E COEFFICIENTS 
C M= I NERT I A COEFFICIENT 

PER=PEP I OD IN SECONDS (CHART PERIOD/CHART SPEED) 

C D=DR AG COEFFICIENT 

REMF=REMA IMDER FUNCTION 

A I » BI =F OUR l E R COEFFICIENTS 

N =NUM3 E D OF DATA SETS 

CN=CGNV ER T I ON FACTOR 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC: 
DIMENS I ON FORCE ( 36 ) 

N = 5 

DO 100 1=1, M 
G=32. 174 

READ ( 5 , 55 ) ( FORCE ( K ) , K=l,36) 

RHG=6 ? . 4/G 
PI =3. 14159 
CNU=0. 0000105 
CL=1 .4946 

INITIATE DATA SETS 



CCCCCCCCCCCCCCC 



r 



READ IN PARAMETERS WHICH CHARACTERIZE THE DATA SET 



R E AD ( 5 , 10 )DIA,AMP , PER , NC ARD , UMX , CN, FMMM , NRUN 
COMPUTE FORCE COEFFICIENTS 

Zl = (2*PER*P=M/(P!*PI*Pi*DIA*D I A * C L * R H 0 * A M P ) 
Z2=(-3-PER*°ER) / ( S * P H 0 * D I 4*P I *C.L * AMP* AM° I 
Z2LS=(-4*PER*?tR ) / ( 3 * P I * P I*RH3*QIA*CL*Ai’lP*AMP) 
C COMPUTE BETA 

BETA=2*PI*AMP/DI A 
R E YNO= ( UMX*D I A ) /CNU 

WRITE (6 , 1 5)0 I A, AMP, PER, UMX, CN,NRUN 
WRI TE ( 6 , 20 ) 

T I ME=0 . 0 
SUM=0. 0 

DO 65 K = 1 , NCARD 
SUM=SUM +FORCE ( K) 

65 CONTINUE 

FMEAN=SUM/NC ARD 
DO 85 K=l, NCARD 
FORCE (K )=FCRCE(.< ) -F ME AN 
85 CONTINUE 

F M MM= F M MM - FM E A N 

CM=0 . 0 

CD=0. 0 

C M L S = 0 

CDLS=0 

F AA= 0.0 

FBB=0 . 0 

FCC=0.0 

FDD = 0 .0 

FEE=0 . 0 

DELTAT=0. 027972020 
DO 200 K = 1 , NCARD 
F=CN*F ORC E ( K ) 

ALPH A = 2 *P I *T I Mr 
S INA = S IN( ALPHA ) 

COSA=CCS( ALPHA ) 

FSINA=OELTATvr*SINA 
F COS a = DC LT at *F*C OS A 

fls=ofltat*f*c:)Sa*abs( cos a ) 

FAAS = 2* PI *F*Dt LTAT*COSA*CDSA*COSA*CQSA*F 
F8BS = 2*P l*F*F< CDSA*ABS< COSA ) * D E L T A T * F 
FCC$ = 2* P i*F*SI NA*COSA*AbS(COSA )*DtLTAT*F 
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FDDS=2*P I*F*SI\'A*SINA*DELTAT*F 

FEES=2*PI*F*F*SINA*DELTAT*F 

F AA=F A AS* FAA 

FBB=FBBS*F88 

F CC= FCCS + FCC 

FDD=FDDS*FDO 

F£E = FE FS * FEE 

C M= F S I M A + C M 

CD=FCOSA+CD 

CMLS=FSIHA *CMLS 

CDLS=FLS*CDLS 

WR I T £ ( S , 30 ) F I ME, ALPHA, EOS A, SIVA , F , FCOSA , FSI NA 
T I ME =T I ME*DELT AT 
200 CONTINUE 
CM = Z 1 * C M 
C0=Z2*C 0 
CMLS= Z 1 *C MLS 
CDLS=Z2LS*CDLS 

CMFF=< Z 1/2 .0) *< FEE*F AA-FCC*FBB ) / ( FDD*FAA-FCC*FCC ) 

CDFF= ( -4.0*Z2/!3.0*PI ) )*( FEE*FCC-FDD*FB3 )/( FDO*FAA-FCC 
*FCC ) 

WRI TE ( fc , 3 5 ) 

WRITE( 6,40)CM,:0,3ETA,REYM3,CMLS ,CDLS , CMFF , CDFF 

ANGLE=0 .0 

WRITE(6,A5) 

T I ME=0 . 0 

DO 300 K= 1 , MCARD 
F=CN-FORCE(K) 

THETA 1= ( ( 2 .O* 3 I ) / 3 60)* ANGLE 

Cl = ( ABS(CCS( THETA 1 ) ) ) *COS ( T HE T A 1 ) 

C 2= RHO* ( { U ‘1 X * * 2 . 0 ) / 2 . 0 ) *DI A*CL 

C3=( (PI*- 2 )OI A*SIN< THETA 1 ) >/(UMX*PER) 

F 1 = ( C.M*C3-CD*C1 ) 

FLS=(CMLS*C3-CDLS*C1 ) 

FF0R=(CMFF*C3-CDFF*C1 ) 

F=F/C2 

REMF = AB S ( F ) -ABS ( FI ) 

FMAX = FMMM*C\'/C2 
REMF=R5MF /FMAX 
RLS=(A8S(F)-ABS( FLS) l/FMAX 
RFOR=( ABS { F) -ABS ( FFOR ) ) /FMAX 

WRITE( 6 , 5 0 ) T I M E » r , c 1 , R E M F , FLS, RLS, FFOR, RFOR 
ANGLE = ANGLE* 10.06993 
TIME=TI ME+DELTAT 
300 CONTINUE 
100 CONTINUE 

10 FORMAT ( 3F9.4,I8,3F8.4,I6) 

1 5 FORMAT ( «1',2X,'0IA = ',F3.4,2X,'AMP = ',F3.A,2X, , PE : ! = », 

1 F8 . A , 2 X , * UMX = « , F8 .4, 2X , • CM- ' , F 8 . 4 , 2 X , ' NRUN = ',16) 

20 FORMAT ( « 0 ' , 3X , • T I ME/ PE R* , 7X , • ALPHA' ,7X,’CCSA',15X,'SIN 
3' ,2X,'F',7X,'F:0SA' , 5 X , 1 F S I N A ' ) 

30 FORMAT ( '0«,3F12.4,Fiq.4,3 c 12.4) 

35 FORMAT (•0',7X,'CM=',7X,'C0=',7X,'6ETA=',7X,'REYNO=',7X 
1 ' C ML S = ' ,7X,'CDLS = ' ,7X, ' C M F F = • , 7X, ' CDFF = • ) 

40 F ORM AT ( • 0 ' , 6 F 1 2 . 4 , 2F 1 2 . 4 ) 

45 FORMAT ( ' 0 • , 7 X , 1 T I M c ' , 9X , ' = ' , 9X , • F 1 ' , 9X , • RE MF • , 9X , ' FL S' 
1 ' RLS ' » 1 2X , ' F FOR ’ , 12X, ' RFOR' ) 

50 FORMAT! 'O' , 6F 1 2 . 4 , 2F 1 5 . 4 ) 

55 FORMAT! FIO. 4 ) . . 

STOP 

END 
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APPENDIX B 



(LIFT COEFFICIENT PROGRAM) 



ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C LIFT COEFFICIENTS 

C H = GR APH ED WAVE HEIGHT (MM) 

C E=ELEVATIGN CALIBRATION (FT/MM) 

C A= AMPL I TUDE (FT) 

C 0=CY L I NOE R DIAMETER (FT) 

C X = CYL I NQER LENGTH (F-T) 

C B = BE T A (DIMENSIONLESS PERIOD PARAMETER) 

C REYMO=R EYNGLD* S NUMBER 

C U'1 = MAX I MUM WAVE VELOCITY (FT/SEC) 

C T=WAVE PERIOD (SEC) 

C FMU= GRAPHED FORCE AWAY FROM BOUNDARY ( MM ) 

C FMD= GRA PHE D FORCE TOWARDS 3 CUM DA RY ( MM } 

C RHO=DENSITY ( SLUGS /CJo I C FT) 

C CLU=L I FT COEFFICIENT AWAY FROM BOUNDARY 

C CLD=L I FT COEFFICIENT TOWARDS BOUNDARY 

C ATTF=FORCE ATTENTUAT I ON COEFFICIENT 

C ATTA = AM PL I TUDE ATTENUATION COEFFICIENT 

C CN^FORCE CALIBRATION (L3/-IM) AT ATTF = I00 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
WRI TE (6 , 200) 

M = 2 

DO 30 K = 1 ,M 

READ (5,100) E,D , X, T , RHO, CM , N» J 
DO 20 1 = 1 ,N 

READ ( 5,400) H, ATTA, FMU,FMD, ATTF 
A = (H/2) *E=M ATT4/ 100) 

UM=2*3. 14 1 7*A/T 
B=UM*T/D 

REYN0=UM*D/0 .00001 05 

CLU=( FMU*ATTE*CN )/ ( . 5 *RHO* J M*UM*0*X* 1 0 0 ) 

C L D = ( F M D * A T T F v C \ )/ ( . 5*RH0*UM* JM *D*X* 1 03 ) 

WRITE (6,300) J , D , B , RE YNO , C LU , C L D 
J = J + 1 

20 CONTINUE 

30 CONTINUE 

100 FORMAT ( 6F8 .A , I 3 , 13 ) 

200 FORMAT (//, 3X ,' R UN ' , a.- X , ' 0 I AM • , 5X , * 5 E T A « , 6X , 1 R E YNO • , 5X 
, « CLU' , 16X,'CL0, / /) 

300 FORMAT* // , 4X , ! 2 , 6X,F6 . A, 3X, F 5. 2, 4X,F8. 1 , 3X, F6.4, 2X,F6. 
400 FORMAT ( F8.4,F8.4, F8.4vF8.4,e8.2) 

STOP 

END 
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